Local density of states of a rf-wave superconductor with inhomogeneous 

antiferromagnetic correlations 
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The tunneling spectrum of an inhomogeneously doped extended Hubbard model is calculated at 
the mean field level. Self-consistent solutions admit both superconducting and antiferromagnetic 
order, which coexist inhomogeneously because of spatial randomness in the doping. The calculations 
find that, as a function of doping, there is a continuous cross over from a disordered "pinned 
smectic" state to a relatively homogeneous d-wave state with pockets of antiferromagnetic order. 
The density of states has a robust d-wave gap, and increasing antiferromagnetic correlations lead to 
a suppression of the coherence peaks. The spectra of isolated nanoscale antiferromagnetic domains 
are studied in detail, and are found to be very different from those of macroscopic antiferromagnets. 
Although no single set of model parameters reproduces all details of the experimental spectrum in 
Bi2Sr2CaCu2 08, many features, notably the collapse of the coherence peaks and the occurence of a 
low-energy shoulder in the local spectrum, occur naturally in these calculations. 



I. INTRODUCTION 

Nanoscale inhomogeneities have been widely ob- 
served in the high temperature superconductor (HTS) 
Bi2Sr2CaCu20s(BSCCO), primarily through scanning 
tunneling microscopy (STM) experiments performed in 
the superconducting stateiS^iiSiSii At present, the ori- 
gins of the inhomogeneity are not well understood, al- 
though it is very plausible that they are directly corre- 
lated with variations in the local doping concentration. 
In BSCCO, the hole concentration is controlled by the ad- 
dition of interstitial oxygen atoms which appear to reside 
w 5 Aabove the conducting CuC>2 layers. Because of the 
short distance, large spatial fluctuations of the Coulomb 
potential are expected in the Cu02 layers, especially in 
underdoped samples where screening is poor. In STM ex- 
periments (which do not measure local doping directly), 
nanoscale inhomogeneities are manifested most strongly 
in the magnitude At of the superconducting gap in the 
tunneling spectrum. Interestingly, regions with small At 
exhibit large coherence peaks at the gap edge, while the 
coherence peaks are essentially missing in regions with 
At > 65 meV. These latter regions are assumed to 
represent an underdoped "pseudogap" phase which may 
be quite distinct from the small-gap "superconducting" 
(SC) regions. Though speculative, this labelling is sup- 
ported by the fact that the "pseudogap" regions occupy 
a large fraction of the strongly underdoped samples, and 
relatively little of the optimally doped samples^ 

The simplest model of the inhomogeneities is that the 
pairing energy depends strongly on doping^SiSiifl with the 
"pseudogap" domains corresponding to regions of large 
pairing energy which are nonsuperconducting because of 
phase fluctuations. It is also possible that secondary 
phases may coexist with the SC state, forming spon- 
taneously, or perhaps being nucleated in hole-poor do- 
mains. Theoretical calculations from the early days of 
high-temperature superconductivitj*ii*i2iiii4 suggest the 
possibility of self-organized stripe formation, and in some 
materials (notably La2- x -j,Nd a Sr x Cu04i^) there is solid 



evidence for stripes, though it has generally been hard to 
substantiate in other materials. A large variety of other 
competing or coexisting phases have been discussed since 
the discovery of HTS, and the list includes charge density 
wave (CDW)fiSiii spin Peierlsr&iS antiferromagnetic or 
spin density wavej22i2ii22i22i pair density waveji^Si, stag- 
gered flu5*2^, and orbital curren t 18 i 26 i 27 i 28 phases. 



Experimentally, optimally-doped BSCC( 
pears to support a fairly straightforward d-wave 



ap- 
BCS 



picture of superconductivity rather we ll f 2ii22i2&2ii2£ al- 
though similar experiments^ have been interpreted in 
terms of commensurate stripe formation with a period- 
icity of « 4ao where ao ~ 5 Ais the lattice constant. 
At lower doping, other recent work£ finds weak "checker- 
board" charge modula tion s with a periodicity close to 
that measured in Refs. [33, with the weight of the mod- 
ulations being reduced as the doping increases. However, 
the situation is not transparent since the modulations 
are only seen at energies larger than the gap edge (con- 
trary to what one might expect in a stripe scenario), are 
only seen in the "pseudogap" regions, and the modulation 
wavelength is comparable to the typical size of the "pseu- 
dogap" domains. There is a further ambiguity in deter- 
mining what, if any, ordering is present: many ordered 
phases (eg. antiferromagnetism) which may be relevant 
to BSCCO couple weakly to the local charge density and 
are not easily identified in STM experiments. Thus, it is 
not clear whether the weak charge modulations seen in 
experiments are the dominant ordering, or whether they 
are secondary manifestations of some hidden order. For 
these reasons, it may be difficult to detect and study co- 
existing order based on spatial modulations of the local 
density of states (LDOS) alone. 

The goal of the present work is to look for signatures of 
inhomogeneously coexisting order in the energy depen- 
dence of the local spectrum. For definiteness, I adopt 
a model in which antiferromagnetic (AF) correlations 
compete with SC order. From a calculational perspec- 
tive, this is the simplest and least ambiguous choice, al- 
though other order parameters — particularly CDW order 
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(including checkerboard order) — are also potentially rel- 
evant to the STM experiments cited above. Many of 
the results of this paper will actually apply broadly to 
other forms of competing order, and I will try distin- 
guish these results from those which are specific to anti- 
ferromagnetism. Having said this, I want to remark that 
antiferromagnetism is a natural choice to make given the 
proximity of the AF and SC phases in the HTS phase 
diagram, and that other authors have studied similar 
modelsii2i2Si2iiS^ In addition, there is mounting evidence 
that glassy (short-ranged) quasistatic AF correlations are 
significant in underdoped HTS; 36 i 37 i 38 i 39 i 40 i 41 including 
BSCCOf 41 " and it is important to understand how these 
correlations are manifested in the LDOS. 

The paper is organized as follows: In Sec. Ill Al I intro- 
duce the model, and perform calculations for a finite- 
sized, inhomogeneously doped d-wave superconductor 
with competing AF and SC order. Short range AF or- 
der arises naturally in the current work because the sys- 
tem is doped inhomogeneously by charged out-of-plane 
donors and AF moments form preferentially in under- 
doped regions. At low doping levels, I find that the 
sclf-consistently determined electronic state resembles 
a "pinned smectic" in which superconductivity is pro- 
nounced along domain walls of the AF background. At 
higher doping, there is a crossover to a fairly homoge- 
neous d-wave SC state with occasional pockets of AF 
order. In all cases, there is a well defined d-wave gap in 
the spectrum. Since the spectral energy resolution suffers 
from finite-size effects, I discuss the LDOS in the context 
of a single underdoped pocket embedded in a homoge- 
neous d-wave superconductor in Sec III Bl Several spec- 
tral features measured in Q , notably the suppression of 
coherence peaks, the appearance of shoulders in the spec- 
trum, and the homogeneity of the low energy spectrum, 
can be understood in these calculations, although no sin- 
gle parameter set reproduces simultaneously all the ex- 
perimentally measured spectral features. One of the most 
important conclusions of this section is that, because of 
the nonlocality of quasiparticles, the local spectrum of 
an AF pocket resembles neither that of macroscopic an- 
tiferromagnets or superconductors (nor is it an average 
of the two): The introduction of inhomogeneity on the 
nanometer length scales leads to a qualitatively new spec- 
trum. This is a significant finding since one of the main 
arguments against coexisting secondary phases is that, 
apart from the special case of a nested Fermi surface, 
any macroscopic ordering which is commensurate with 
the lattice has a spectrum which is not particle-hole sym- 
metric, in contradiction with experiments. I find, how- 
ever, that inhomogeneous ordering on nanometer length 
scales may, in fact, yield a particle- hole symmetric spec- 
trum. These calculations are interpreted in terms of a 
three-band model of homogenously coexisting SC and AF 
order in Sec. Ill CI The issue of how charge modulations 
arise in this model is discussed in Sec. Ill Dl Conclusions 
are presented in Sec. IIIII 



II. CALCULATIONS AND RESULTS 

A. Inhomogeneously doped superconductor 

The basic Hamiltonian is the Hubbard model with a 
long-range Coulomb interaction and SC pairing interac- 
tion: 

i ijij 
ij 

where Cj a is the spin-rr annihilation operator at site j, hi a 
and hi are the spin-resolved and total charge density op- 
erators at site i, and r, is the position of the i th site. I use 
a third-nearest neighbor conduction band with param- 
eters to, ..A3 describing the on-site potential, nearest, 
next-nearest, and third-nearest neighbor hopping ampli- 
tudes. Throughout this work, all energies are given in 
units of \ti\ which (for reference) is ~ 0(100) meV. I 
take {ti,t 2 ,t 3 } = { — 1,0.25,-0.1} and adjust to to give 
the desired filling. The long range Coulomb interaction is 
V(r) = (e 2 /eao)|r| -1 where r is measured in units of the 
lattice constant do, e 2 /eao = 1, and the on-site interac- 
tion is absorbed into the Hubbard U term: V(0) = U/2. 
The impurities are located at positions R, which sit a 
distance d z = 1.5ao above randomly chosen lattice sites. 
The final term in the Hamiltonian is added as an ansatz 
to describe SC order arising from spin-interactions be- 
tween neighboring sites. The local bond order parameter 
Aij — —j}{cj[Ci-\ + cj^Cii) i s determined self-consistently 
for nearest neighbor sites i and j. The Coulomb in- 
teraction is treated in the Hartree approximation, and 
the effective mean-field Hamiltonian can be diagonalized 
numerically to extract eigenstate wavefunctions and the 
corresponding eigenenergies. The fields Ay and Tiifj are 
iterated to self-consistency on small lattices with between 
20 x 20 and 40 x 40 sites. The calculations are uncon- 
strained, and are seeded with a finite antiferromagnetic 
moment. In order to improve convergence, which is prob- 
lematic when magnetic moments form, a combination of 
Thomas-Fermi and Pulay method charge-mixing is used 
at each iterative stepi 4 ^ 

Figure ^ shows typical results for an underdoped su- 
perconductor for the self-consistent hole density um = 
1 — m with Hi — (hi), the staggered AF moment rrii = 
ig»Q ri^ n ^ — Tin) w ith Q = ( 7r J 7r )i an d the d-wave order 

parameter A^ = \ l) 5 " \i+s where S is summed 

over the four nearest-neighbor sites. The average hole 
density is rih = 1 — n « 0.07 holes/site, but the hole dis- 
tribution is quite inhomogeneous. In a homogeneously 
doped sample with J — 1.5 and U — 3.2 (these are typ- 
ical for this work), there is a first order phase transi- 
tion between SC and AF phases at a hole doping level 
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FIG. 1: Self-consistent solutions of the model, (a) Charge 
density, (b) staggered magnetization, and (c) d-wave gap are 
shown for a 30 x 30 lattice with 35 donor impurities of charge 
Z = — 2e. The model parameters are U = 3.2 and J = 
1.5. The corresponding hole doping level is nn ~ 0.07. The 
staggered magnetization for a single underdoped disk of radius 
4ao is shown in (d). Notice that a small incommensurate 
moment is induced outside the underdoped disk. The LDOS 
along the line from R to B is shown in Fig. [3] 



Uh = 0.07. In the inhomogeneous system, the situation is 
more complicated. The spin polarization saturates near 
its bulk value in undoped regions whose diameter ex- 
ceeds £af ~ ti/Um, where m is the staggered moment. 
In smaller underdoped regions, the staggered moment is 
roughly proportional to the diameter of the region. It is 
worth stressing that this behavior is very different from 
single-phase models in which the magnitude of the lo- 
cal order parameter is directly correlated with the lo- 
cal charge density^ regardless of the size of the domain. 
Note also that, although the doped and undoped regions 
in Fig. n lie firmly on either side of the first order phase 
transition separating AF and SC phases, both order pa- 
rameters are finite throughout the system because of a 
pronounced proximity effect. In this sense, the intro- 
duction of disorder in the doping leads to a qualitative 
change in the phase diagram. This aspect of the calcu- 
lations appears to be consistent with neutron scattering 
studies in LSCC 4Q suggesting that AF and SC coexist lo- 
cally. One factor which appears inconsistent with exper- 
iment is that the d-wave order parameter is suppressed 
by static AF correlations, whereas there is good evidence 
that it actually grows rapidly as the insulating phase is 
approached in HTS. This disparity may be the result of 
the simplicity of the current mean-field approximation. 
In a more sophisticated treatment, the suppression of SC 
order will be compensated to some extent by the fact that 
the pairing interaction is doping dependent: Numerical 
studies of the t-J mode&ia find that J ~ (1 + n h )^ 2 . 
Since this result was originally derived for homogeneous 




FIG. 2: Doping dependence of the AF and SC phases. The 
three rows display the hole density (top), SC gap (middle) 
and average density of states (bottom) for three different 
donor-impurity concentrations. The columns correspond to 
20 donor-atoms (left) 35 donor-atoms (middle) and 70 donor- 
atoms (right). Contours show the domain walls of the stag- 
gered AF moments. Parameters are as in Fig. □ The color 
scales are identical for all panels within a row, and are the 
same as in Fig.^a) and (c). 



systems, and cannot be trivially extended to inhomoge- 
neous systems (but should not change our conclusions 
qualitatively), I will defer its discussion rather than in- 
troduce an ad hoc local renormalisation of J. 

Figure [21 shows how the coexisting phases evolve with 
hole doping. At low doping, the situation qualitatively 
resembles a pinned smectic. Smectic phases have been 
proposed as a natural mechanism by which doped antifer- 
romagnets can accomodate holes while minimizing both 
the hole kinetic, antiferromagnetic exchange, and long 
range Coulomb energies42, Consistent with the smectic 
picture, the AF moments in Fig. [21 spontaneously form 
7r-shifted domains whose boundaries are pinned to donor- 
impurity locations, and the SC order parameter is 
largest along the domain walls. However, there are a 
number of differences between the current "weak cou- 
pling" mean-field calculations and the canonical "strong- 
coupling" smectic picture. First, because of frustration 
introduced by next-nearest neighbor hopping, the AF 
phase is never fully polarized. This same frustration leads 
to a gapless quasiparticle spectrum in the pure AF phase 
provided U is less than some model-dependent critical 
value. As a consequence, holes are not confined to do- 
main walls but are mobile throughout the volume of the 
sample. In other words, although there are static AF cor- 
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relations, the system is on the metallic side of the metal- 
insulator transition. A further consequence, which dis- 
tinguishes weak and strong-coupling approaches, is that 
the SC order parameter remains finite everywhere in the 
weak-coupling calculations. 

As hole doping increases, the AF phase is suppressed 
through a proliferation of domain walls. In Fig. [21 signif- 
icant AF moments form at higher doping only in regions 
where, due to randomness in the donor-atom distribu- 
tion, undoped regions have diameters larger than £af- 
The system is then better understood as consisting of 
isolated AF pockets embedded in a relatively homoge- 
neous d-wave superconductor. In this model, annealing 
(which tends to homogenize the charge distribution) will 
have significant effect on the extent of AF order. 

Figure [21 also shows the spatially averaged density of 
states p(ui) for each of the disorder configurations. Al- 
though the data is noisy, several clear features are evi- 
dent: First, as one underdopes, there is a gradual sup- 
pression of spectral weight on an energy scale which is 
large relative to the SC gap. Second, there is a robust 
rf-wave gap, even in situations where a large fraction of 
the sample is antiferromagnetic. I emphasize that this 
latter effect is not necessarily anticipated since, in the 
absence of a nested Fermi surface, AF ordering tends 
to destroy the particle-hole symmetry of the spectrum. 
Third, as one underdopes, the superconducting coherence 
peaks are suppressed. As discussed in the introduction, 
the suppression of coherence peaks is one of the hallmarks 
of the pseudogap phase of the underdoped cuprates. To 
my knowledge, this is the first reproduction of such an 
effect in terms of a static mean field model. 

These results are rather encouraging and, ideally, the 
next step should be a detailed examination of the LDOS. 
However, finite size effects limit the spectral resolution 
to the extent that the LDOS is impossible to interpret. 
Consequently, I will focus for the remainder of the pa- 
per on intermediate doping levels where one can improve 
the energy resolution by studying isolated underdoped 
pockets embedded in a large superconducting domain. A 
more detailed exploration of the pinned smectic phase re- 
quires a different approach and is, unfortunately, beyond 
the scope of this work. 



B. Single underdoped pocket 

It is difficult to discuss the STM spectrum in detail 
for finite-sized lattices because of the discreteness of the 
spectrum. For a 30 x 30 lattice, one might typically have 
~ 100 subgap states, with the resulting spectrum being 
too noisy for anything other than the grossest analysis. 
I therefore study a single, isolated, AF pocket which is 
embedded in a homogeneous background potential cor- 
responding to a hole doping level of p « 0.15 (in fact, 
the hole doping level is less important than the Fermi 
surface shape, and should not be taken too seriously). 
In this calculation, a positively charged disk of radius R 
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FIG. 3: Local density of states, (a)-(d) The LDOS is shown 
for different model parameters. The spectra, offset for clarity, 
are taken at a sequence of sites extending radially outwards 
along the (010) direction from the center of an isolated un- 
derdoped disk of radius R. The path along which spectra are 
measured is shown in Fig. H]d), with points B and R corre- 
sponding to the blue (top) and red (bottom) curves in this 
figure. The heavy black curve in each panel indicates the site 
at which the staggered magnetization falls to half the max- 
imum value m max . Diamond symbols indicate ±2A' d ' (the 
estimated coherence peak energies) at each site. The first 
three panels are for (a) U = 3.2, R — 1.5ao, (b) U = 3.2, 
R = 4.0a , and (c) U = 3.4, R = 4.0a . J = 1.5 throughout. 
In (d), spectra are calculated for a non-self-consistent model 
of a pure antiferromagnetic pocket of radius R = 6.0 with 
U = 3.4 embedded in a pure d-wave superconductor. The 
heavy black curve marks the sharp boundary between the 
AF and SC domains. For comparison, the DOS for homo- 
geneously coexisting AF and SC order are shown in (e)-(h). 
Cases are (e) (M, A (d) ) = (0,0.3), (f) (0.3,0.3), (g) (0.6,0.3), 
(h) (1.0,0). 



sits d z = 1.5ao above the conducting layer. The charge 
on the disk is adjusted so that the site under the center 
is half-filled. The charge, magnetization and SC gap are 
calculated self-consistently, and an example of the self- 
consistent magnetization for a disk of radius R = 4ao is 
shown in Fig. ^d) . In order to obtain a high spectral 
resolution, the underdoped pocket is embedded in a ho- 
mogeneous 200 x 200 region and a recursion technique 44 
is used to calculate p(r,w). In this way, I avoid spurious 
structures associated with the discreteness of the spec- 
trum on finite lattices. Figure [31 shows the LDOS along 
cuts through the centre of an AF pocket for different val- 
ues of {U, R}. 

The spectra in Fig.[31show a smooth evolution from re- 
gions where the ci-wave order parameter is dominant, to 
the central region where AF correlations are large. The 
particle-hole asymmetry at large energies comes from a 
van Hove singularity at lu sa — 0.5|ii|. At lower energies, 
the spectrum is determined by the interplay between SC 
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and AF order. There are several noteworthy features of 
this calculation. First, there is an overall suppression of 
spectral weight at site i on an energy scale M, = Urrii. 
For a nested Fermi surface, with homogeneous magneti- 
zation, Mi is the energy of the AF gap, but in the ab- 
sence of nesting, AF correlations lead to a shift of states 
away from the Fermi level even when a true gap does not 
open. This shift is a precursor to the formation of lower 
and upper Hubbard bands, and is consistent with the ob- 
servation of spectral weight shifts on a large energy scale 
as a function of doping in the cuprates. Two factors tend 
to suppress the nucleated moment: the frustration intro- 
duced by the next-nearest neighbor hopping (ie. the ab- 
sence of complete nesting of the Fermi surfaces) , and the 
competition with the superconducting phase. I want to 
emphasize a consequence of this which may not be intu- 
itive: Although a large U may not generate a substantial 
moment to, the energy scale M over which the quasiparti- 
cle spectrum is affected by magnetic correlations can still 
be quite large. This is likely to be a universal feature of 
models of competing order in the cuprates. One could 
similarly imagine that, in a model with a charge ordered 
phase, frustration due to imperfect Fermi surface nesting 
and competition with SC order would tend to suppress 
the magnitude of charge modulations but still affect the 
spectrum over a relatively large energy scale. 

A second feature of Fig. is that the local dispersion 
in the AF pocket near the Fermi level is quasilinear when 
Mi < A| d) [Fig. EJa)]. This is a significant result since 
one of the main arguments against coexisting commen- 
surate order in the cuprates is the absence of a linear dis- 
persion at the Fermi level. When M t > A (rf) , the LDOS 
becomes particle-hole asymmetric: as Mi increases (ei- 
ther as one moves into the AF pocket, or as one turns 
up U), a shoulder develops in the dispersion at low en- 
ergies, which ultimately evolves into a well-defined reso- 
nance [Fig. HJb)]. The energy of the resonance depends 
on the details of the band structure, and on A"', but is 
a universal feature in nearly all numerical results. 

One of the most interesting aspects of Fig. is the 
evolution of the coherence peaks between the SC and AF 
regions. There are actually two qualitatively different 
ways in which this occurs. In Figs. |^ a) and (b), the co- 
herence peaks shift to lower energies, sharpen, and lose 
spectral weight as one moves into the AF domain. The 
coherence peak positions approximately reflect the local 
value of A^, which is what one might naively expect for 
a smoothly varying Hamiltonian. The situation is differ- 
ent in Fig. GUc) where the magnetization is larger: the 
coherence peaks (starting at a point exterior to the AF 
domain and moving inwards) rapidly collapse, but shift 
rather little. The absence of a shift indicates that one is 
seeing the decaying tails of bulk BCS-like states. In other 
words, antinodal quasiparticles from the SC domain tun- 
nel (rather than propagate freely) into the AF domain, 
and decay over some characteristic distance which de- 
termines the extent of the coherence peaks into the AF 
domain. I will argue below that this arises from a mis- 



match in the SC and AF energy dispersions: when M is 
sufficiently large, the states at the antinodal /c-vector are 
gapped in the AF domain. 

The fact that the coherence peaks sharpen as one 
moves into the AF region in Fig. I^a) and (b) is the 
result of the fact that the Fermi surfaces nest at isolated 
points (the contours ek = and ek+Q = shown in 
Fig. 01 intersect at two points). Because of this pecu- 
liar nesting, spectral weight is removed at energies both 
above and below the antinodal saddle point energy (the 
point which generates the coherence peaks) by the AF 
correlations, but the saddle point itself survives until the 
AF moment becomes very large. Thus, the coherence 
peaks lose weight by narrowing rather than by being sup- 
pressed. Other competing phases (such as charge density 
waves) which nest differently should have a qualitatively 
different effect on the coherence peaks. 

One surprising aspect of Fig. EJb) and (c) is that al- 
though the transition from SC to AF domains occurs dif- 
ferently depending on the magnetization, the spectrum at 
the core of the AF pocket is quite similar. For compari- 
son, a non-self-consistent calculation is shown in Fig.^d) 
for the ansatz 

_ f 0, |r i |< J Ror|r i |< J R 
13 \ 0.3, otherwise 

to = f 0.3, |r 2 | < R 
1 \ 0, otherwise ' 

with R = 6a,Q. Again, the spectrum at the core of the 
AF disk is quite similar to that of the self-consistent cal- 
culations shown in Fig. Etb) and (c) . (note the similar- 
ity in the peak positions), but bears little resemblance 
to the spectrum of the macroscopic AF phase shown in 
Fig-Eth). This calculation demonstrates that the bound- 
ary conditions (the coupling between the AF pocket and 
the SC bulk) have as large an impact on the local spec- 
trum at the core of the AF pocket as the local value of 
the SC order parameter. This is a central result: when 
the scale of the inhomogeneity is atomic, one cannot as- 
sume a direct correspondence between the local ordering 
and the local spectrum. In the following section, I argue 
that a qualitative understanding of the inhomogeneously 
doped system can be developed from a model of homo- 
geneously coexisting SC and AF order. 



C. Homogeneously coexisting order 

I consider a three-band model of a homogeneous sys- 
tem with coexisting SC and AF long range order. The SC 
order parameter has the usual form Ak = A^ d '(cosk x — 
cosfcy) and the AF order parameter is M with to.; = 
M/U. I adopt the same dispersion as before, with = 
t + 2ii(cos k x + cos k y ) + 4t 2 cos k x cos k y + 2i 3 (cos 2k x + 
cos2k y ) and to — 0.7|ti|. For a complex frequency 
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FIG. 4: Constant energy contours of the coexisting AF/SC 
model. Parameters are the same as Fig.^c). A few contours 
are shown for positive energies for the upper (dashed red) and 
middle (solid blue) bands. The zero-energy contours of ek and 
£k+ q are shown for reference. For small energies, the contours 
have the "banana" -like shape expected for ci- wave supercon- 
ductors, as well as an antiferromagnetic shadow band. There 
is a saddle-point singularity at uj = 0.38 (labelled "A" in the 
figure) which marks the end of the linear dispersion, and gives 
rise to the ui = 0.38 van Hove singularity in Fig.^f). There is 
a second saddle-point at B in the figure, which gives rise to the 
superconducting coherence peaks. The energies of the saddle- 
points depend on both A^ and M. For sufficiently large M, 
the saddle-point at B evolves into a simple band minimum, 
corresponding to the lower edge of the upper Hubbard band. 



z = UJ + i0 + , the Green's functions satisfy 
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where k = k Q, and where G kk ,(uj) and F kk ,(ui) 
are Fourier transforms of the retarded and anomalous 
Green's functions: 

Gkk'W = -i(Mt),ci, a (0)}}@(t) 

Fkk'W = -*<{cl k ^),4v(o)}>e(t). 

The density of states p(w) — — -Im J2k g G kk (uj + i0 + ), 
plotted in Fig. Ete)-(h), is determined by the poles of 
G kk (cj). For reference, a few constant-energy contours 
of the spectrum are shown in Fig. 0] for a case with 
M = A^ = 0.3. The interested reader is directed to 
Ref. for an extensive discussion of the normal state 
spectrum of this model. In the coexisting state, there 



are two features of interest: first, at low energies the 
spectrum resembles that of the pure superconductor and, 
second, there is a new saddle-point singularity (at "A" 
in Fig. which arises because of the coexisting order. 
In the limit M S> A, the origin of the saddle-point is 
fairly transparent: G£ k (z) has three poles correspond- 
ing to upper and lower magnetic bands with dispersion 
E± = (e k + e k )/2 ± [(e k - e^ 2 /4 + M 2 ] 1/2 and a hole- 
band with dispersion E = —e k . When Ak is nonzero, 
there is an avoided crossing of E and E_ which results 
in the saddle point. Both features of the dispersion are 
evident in p(u>) [Fig. EJf)], which resembles the pure d- 
wave superconductor at low uj, and has a resonance at 
the saddle-point energy, uj = 0.38. This model appears 
to capture several aspects of the inhomogeneous spectra 
in Fig. EJa)- (d). First, it predicts the overall suppression 
of spectral weight on magnetic energy scales. Second, it 
predicts the occurence of a subgap resonance. Third, it 
shows that while a d-wave-like tunneling gap survives to 
fairly large values of M, there is an inward shift in the po- 
sition of the apparent "coherence peaks" as M increases. 
This is purely a band structure effect resulting from the 
reduction of spectral weight at the antinodal points, and 
has nothing to do with a reduction of A^ d \ Fourth, it 
suggests that because AF nesting does little to disrupt 
the band structure near the nodal points, there will be a 
much smaller Fermi surface discontinuity for nodal quasi- 
particles crossing between SC and AF domains than for 
antinodal quasiparticles. This mechanism is one possi- 
ble explanation for relative uniformity of the low energy 
spectrum measured in STM experiments when compared 
with the spectrum near the gap edge. 

This model makes one other, somewhat subtle, predic- 
tion which appears to be relevant to the numerical work. 
Fig-Glshows several constant energy contours for the mid- 
dle and upper bands. For small M, the bottom of the 
upper band lies below the top of the middle band. As M 
is increased, however, a gap will appear in p(u>). Because 
the nesting points (the points at which e k = ek+cj) lie 
near the antinodal points, the gap, when it opens, does 
so near the coherence peak energy, as in Fig.|3Jg). [No- 
tice that the energy at which the gap appears depends on 
both the band structure, and on A^ d ^. Hence, the gap in 
Fig. |HIg) appears at a higher energy than for a pure an- 
tiferromagnet. Furthermore, tends to enhance the 
magnitude of the gap.] The key difference between Fig.|3 
(b) and (c), where the coherence peaks evolve smoothly 
in the former and collapse in the latter, appears to be 
the presence of an AF gap in the spectrum. I remark 
that while there is a threshold value of M at which a gap 
forms, other kinds of order may not have such a thresh- 
old. For example, a charge density wave which nests 
between parallel antinodal sections of the Fermi surface 
may gap out the coherence peaks for any degree of or- 
dering. 

A second consequence of having a gap in p{ui) is that 
scattering resonances may produce exponentially local- 
ized bound states at energies within the energy gap. The 
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FIG. 5: Response kernel for disorder. The imaginary part of 
As(q, uj) is shown for three values of lu for the coexisting order 
(left column) and pure superconducting (right column) cases. 
The parameters are the same as in Fig.^f) (left column) and 
(e) (right column). At low energy, the kernel is similar for 
both models (top row). The saddle point (labelled "A" in 
Fig. causes the peak at q ~ (tv/2,0) to split into two at 
uj = 0.35. 



resonance at w = 0.35 in Fig. Etc), for example, is very 
sharp, and localized to only the few sites nearest the core 
of the AF domain. Furthermore, it's energy is close to 
where the spectral gap opens in the pure AF, making 
it a good candidate for the kind of local resonance dis- 
cussed here. There are several sources of scattering — the 
inhomogeneity of the SC and AF order parameters, and 
the impurity Coulomb potential — which are not included 
in the homogenenous model, which could give rise to a 
bound state. 



D. Weak charge modulations 

Weak charge modulations have been observed in un- 
derdoped BSCCOj£i2S sparking an ongoing debate as to 
the extent to which spatial modulations of the LDOS 
can distinguish Friedel oscillations of quasiparticles from 
a tendency towards charge-ordering. In this section I 
will address two slightly different questions which arise 
from this debate. First, I will discuss charge ordering 
in the context of the current calculations. Second, I will 



discuss the broader question of the extent to which a hid- 
den order (an order which does not couple directly to the 
charge) can be revealed by quasiparticle scattering. 

In Sec. Ill Al the self-consistent calculations show that 
a large charge inhomogeneity occurs in the hole-doped 
domains, and arises because of randomness in the impu- 
rity locations. In contrast, the undoped domains are re- 
markably homogeneous simply because they are free from 
impurities. The charge inhomogeneities are not particu- 
larly evident at energies near the Fermi level, but become 
apparent at large energies. The calculations do not find 
a local charge ordering, although the kind of weak mod- 
ulations seen in Ref. [7j would be difficult to see on the 
finite-size lattices used in numerical calculations. Even 
if the current calculations do not admit charge ordering, 
it is likely that only minor modifications to Eq. Q are 
needed to generate charge-ordered phases. By analogy 
with the self-consistent solutions for the AF phase, one 
would expect that these charge ordered phases would co- 
exist with superconductivity throughout the system be- 
cause of the proximity effect. This is in contrast to ex- 
periments in BSCCOfi however, which see ordered charge 
modulations only in "pseudogap" domains. A resolution 
to this puzzle may be the fact that the charge distribu- 
tion is uniform in underdoped domains, but is disordered 
in hole-doped domains. If the CDW order parameter is 
easily pinned by donor atoms, then it may also be locally 
suppressed by donor-related disorder. 

The second question is whether hidden order can be re- 
vealed through the Fourier transformed density of states 
of a disordered superconductor. In order to get some 
sense of how disorder affects the LDOS, I return to the 
model of homogeneously coexisting order described in the 
previous section. I calculate the response kernelj22i2i 



Aa(q )W ) = ^Tr CT [G£ k (o,)G£ 



k+qk+q 



-^kkM^k+qk+qML 

which describes the effects of scattering from impurities 
on the Fourier transformed density of states p(q, uj). The 
results are shown in Fig. At low ui, Ag(q, uj) is similar 
for both coexisting order and for the pure d-wave super- 
conductor while, at higher energies, the effects of anti- 
ferromagnetism become significant. At u) = 0.35 (the 
energy of the saddle point marked "A" in Fig. the 
pronounced resonances along the (n, 0) and (0, n) direc- 
tions are split by the antiferromagnetism. One surprising 
result of these calculations is that A 3 (q, uj) differs signifi- 
cantly from the pure d-wave result even for energies much 
larger than M, suggesting that the effects of even weak 
ordering should be easily visible. At the same time, how- 
ever, there is no obvious signature of the AF Q-vector 
in the response kernel at most energies. In other words, 
AF ordering distorts the response kernel from the bare 
kernel, but does so in a nontrivial way. In particular, all 
features in As(q, uj) disperse with uj. Thus, it appears 
that unless the nucleated order couples directly to the 
charge density, it will generally be difficult to distinguish 
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different kinds of order from the Fourier transformed den- 
sity of states. 

III. CONCLUSIONS 

I have studied a mean-field extended Hubbard model in 
which charge is doped inhomogeneously because of ran- 
domness in the donor-atom positions. AF and SC order 
compete, and in the homogeneous case, are separated by 
a first order phase transition. Self-consistent calculations 
find that, because of inhomogeneity in the local doping, 
AF order coexists inhomogeneously with superconduc- 
tivity. The AF moments spontaneously form Tr-shifted 
domain walls which are pinned to donor-atom sites. At 
low doping, the self-consistent solutions resemble pinned 
smectics, ie. quasi-one-dimensional superconductors run- 
ning along AF domain walls. Because of the proximity 
effect, both SC and AF correlations are actually present 
throughout the lattice. This picture appears to be con- 
sistent with neutron scattering studies in LSC04& sug- 
gesting that AF and SC order coexist, and naturally ex- 
plains the field-dependence of the AF moment, since any 
suppression of superconductivity by a magnetic field will 
enhance the AF moment. 

At higher doping, the self-consistent calculations 
evolve towards a homogeneous d-wave superconductor in- 
terspersed with underdoped pockets with large AF mo- 
ments. This latter "phase-separated" system superfi- 
cially resembles the situation in BSCCO, although it is 
generally unclear whether the "pseudogap" domains in 
BSCCO have any kind of nucleated secondary phase. 
To address this question, I studied the local spectrum 
of a single, isolated AF pocket embedded in a homoge- 
neous SC background. While no single calculated spec- 
trum reproduces all details of the experimental measure- 
ments, several features such as the collapse of the coher- 
ence peaks, the occurence of low-energy spectral features, 
and the relative homogeneity of the low energy spec- 
trum, are broadly consistent with the kinds of spectra 
measured in, for example, Ref. [3. Certain experimen- 
tal aspects — notably the presence of weak nondispersing 
charge modulations — are not reproduced in my calcula- 
tions. In general, the calculated spectra at low energies 
show a richer spectrum of peaks than is observed exper- 
imentally. 

At this point, the effect of disorder on the spectrum of 
the isolated AF pocket is not understood. Earlier studies 
of point-like defects in <i-wave superconductors show that 
this may not be a trivial effect. A single strong-scattering 



point-like impurity introduces a sharp resonance near the 
Fermi level. As the disorder level increases, the reso- 
nances split, are inhomogeneously broadened, and evolve 
into an impurity band (see Ref. [4r| for a recent sum- 
mary). When the response of the SC order parameter 
to the disorder is included self-consistently, the SC gap 
tends to restore itself by shifting spectral weight away 
from the Fermi level. Indeed, it is a general feature of 
interacting electrons in disordered media that the sys- 
tem can lower its energy by suppressing the density of 
states at the Fermi level. Further calculations, currently 
in progress, are needed to establish whether all the spec- 
tral features discussed in Sec. Ill Bl survive in the disor- 
dered limit. 

Finally, although the calculations were performed for 
a model in which superconductivity and antiferromag- 
netism compete, I expect many of the findings to apply 
to other models of competing order. Three results, in 
particular, are expected to be general. First, when the 
domain sizes are small (as they appear to be in BSCCO), 
the proximity effect is extremely important, and has a 
significant impact on the local density of states. It was 
never found, in calculations, that the spectrum of the 
antiferromagnetic pocket resembles that of a bulk anti- 
ferromagnet. Rather, a better toy model appears to be 
one of coexisting homogeneous superconductivity and an- 
tiferromagnetism. Second, the gapping of the spectrum 
near the antinodal points by local ordering is a mean- 
field mechanism by which coherence peaks may be lo- 
cally suppressed. Up to now, it has been generally under- 
stood that suppression of coherence peaks occurs through 
strong inelastic scattering at higher energies. For AF or- 
der, there is a threshold value of the magnetization for the 
gapping of the antinodal quasiparticles, but this may not 
be universal and other kinds of order (eg. CDW) may lead 
to suppression of coherence peaks for even small order- 
ing. Finally, these calculations also suggest a natural rea- 
son that nodal quasiparticles should be less affected than 
antinodal quasiparticles by charge inhomogneities, since 
the Fermi surface mismatch between domains is smallest 
for the nodal quasiparticles. 
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